h = [86 163 231 288 339 381 413 437 453 458 456 440 400 356];
plot(h,'+')
refcurve([-4.9 100 0])
